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Abstract. In this work we theoretically analyze a circuit QED design where 
propagating quantum microwaves interact with a single artificial atom, a single Cooper 
pair box. In particular, we derive a master equation in the so-called transmon regime, 
including coherent drives. Inspired by recent experiments, we then apply the master 
equation to describe the dynamics in both a two-level and a three-level approximation 
of the atom. In the two-level case, we also discuss how to measure photon antibunching 
in the reflected field and how it is affected by finite temperature and finite detection 
bandwidth. 
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1. Introduction 

In recent years, the field of circuit quantum electrodynamics [TJ|2] (circuit QED) has 
become one of the most promising platforms in the study of light-matter interaction. 
One of the most important breakthroughs in this field was the achievement of strong 
coupling between light and matter, or microwave photons and Josephson-based artificial 
atoms [3111]. Since then, many experiments have been carried out within the framework 
of superconducting circuits [SHS] , revealing a wide variety of novel quantum phenomena. 
Most of these experiments share a common feature, namely the interaction between 
artificial atoms and isolated modes of the electromagnetic field in a cavity. Within circuit 
QED, there is now a growing interest in studying propagating fields interacting with 
artificial atoms, due to e.g. its potential interest in Condensed Matter \X0\ and all-optical 
quantum information [TT]. Theoretically, coherent coupling between a superconducting 
qubit and a one dimensional continuum of modes was discussed some time ago [T2][T3] , 
and there exists now a growing number of experiments investigating this system 
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In this manuscript we report on an in-depth microscopic description of the 
coherent coupling between a field propagating through an open transmission line and 
a superconducting artificial atom based on the single Cooper-pair box (SCB) [22H28] . 
In more detail, we analyze the so-called transmon regime [Hl[29] and study the photon 
transport properties of this system according to different approximations. On the one 
hand, in the two-level approximation, and under certain conditions, the qubit behaves 
as a saturable mirror [T2|, [T3]. On the other hand, including a second excited state of 
the transmon, we can effectively make the medium transparent for the incident photons 
using a coherent control field on resonance with this second transition. Finally, we also 
discuss how the the photon antibunching observed in the reflected field is reduced by 
finite temperature and finite detection bandwidth. Our theoretical predictions are in 
full agreement with recent experiments [T8 |[T9ll2~T] . 

This paper is organized as follows: In section El we derive the master equation 
of a SCB coupled to an open transmission line (TL). In section |2~T| we start from a 
discretized lumped-element description of the TL and in section 12.21 we proceed to the 
continuum description. In section I2.3[ we discuss the regime where the system can be 
described as a SCB weakly coupled to the the voltage of the TL, at the coupling point. 
We thus arrive at the Hamiltonian of a voltage biased SCB, weakly coupled to a bath 
of harmonic oscillators, i.e. the electromagnetic modes of the TL. Making standard 
weak coupling approximations, we then derive a master equation on Lindblad form in 
section 12.41 and attach a coherent drive in section 12.51 For simplicity, we go through 
these derivations considering a SCB at the end of a semi-infinite TL, but in section [2761 
we discuss how the master equation can be straightforwardly extended to an arbitrary 
number of semi-inifinte TLs, all meeting at the SCB. In particular, this includes the 
important case of a single infinite TL. 

In section [3j we then apply the master equation to a few experimentally relevant 
cases [lElIlH] . Section 13.11 is devoted to the reflection and transmission of a single near 
resonant coherent drive, while section [3T21 includes two coherent drives, where one is used 
to control the transmission of the other. Finally in section 13.31 we investigate how the 
photon antibunching observed in the reflected field is influenced by finite temperature 
and finite detection bandwidth. 

2. Model 

In this section we present a general formalism of the light-matter scattering in a one- 
dimensional continuum from a microscopic point of view. We start from a Hamiltonian 
description, arriving at the well-known input-output relations for the microwave field. 
We then follow the usual approach [30] to describe the joint state of the light-matter 
system by introducing dissipation, resulting in the standard quantum optical master 
equation. 
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Figure 1: Discretized circuit describing the interaction of a single Cooper-pair box 
with microwaves photons propagating in a semi-infinite transmission 
line. 



2.1. Discrete Circuit Model, Hamiltonian and Equations of motion 

Consider a semi-infinite transmission line with characteristic inductance Lq and 
capacitance Co per unit length. We discretize the transmission line (31] in units of the 
small length Ax, which we take to zero at the end of the calculation. The transmission 
line nodes are numbered with negative integers, while the SCB island node has index J 
and its Josephson junction has a capacitance Cj to ground and a Josephson energy Ej. 
The SCB is coupled to the transmission line at the zeroth node, through the capacitance 
C c , as depicted in Figured! 

To describe the circuit dynamics we use the node fluxes $ Q (t) = /* dt'V a {t') as 
coordinates [32]. They are the time integrals of the node voltages and although being 
less intuitive than the voltages, this choice greatly simplifies the description of the 
Josephson junction. Starting from a circuit Lagrangian we can derive the discrete circuit 
Hamiltonian [33] 



2Cj 2C C V H J Ax ^ 2C 2L 

J c v 7 n<0 u 

where the charges p a are the conjugate momenta to the node fluxes $ a , fulfilling the 
canonical commutation relations 

where denotes Kronecker's delta. From the Hamiltonian we get Heisenberg's 
equations of motion for the transmission line operators (n < 0) 

a a P« a $n-i - 2$ n + $ n+ i 

= Axcy dtPn = axT ' (2) 

and for the SCB operators 

o~ Po + Pj . Po C s pj 

<9t$o = 7^ — + 7^ = TrprPo + 7^, (4) 
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h \ h 

i 

Cj 



d t $j = , (6) 



where Cy, = C c + Cj. 
2.2. Continuum limit 

In the continuum limit As — > 0, the charge of each transmission line node will 
go zero together with the node capacitance. Thus we define a charge density field 
p(x n ,t) = p n (t)/Ax and a flux field $(x„,t) = where we define the spatial 

coordinate x n = nAx for n < 0, along the transmission line. The continuum equations 
of motion for the transmission line (x < 0) are 

d tP (x) = %p^, d t *(x)=?&. (7) 

These are the equations of motion for the massless Klein-Gordon field, having freely 
propagating left and right moving solutions with velocity v = 1/a/L Co. Therefore, we 
can write the general solution for x < as a linear combination of right and left moving 
second-quantized fields, 



(8) 

where k u — u/v and Zq = ^JLq/Cq is the characteristic impedance of the transmission 
line. The operators aj' annihilate a left /right-moving photon with frequency u, and obey 
the bosonic canonical commutation relations, [a£~, (a,*/)*] = [a£, = 5{oj — u/) and 

[a^,(a^,y] = [a^,a^i] = 0. Finally, we note that in the continuum limit ([3]) changes 
into 

dtPo = F • (9) 

To describe the system dynamics, we first need to specify the incoming, right moving 

field 

$™(t) = $_>((T,t). 

Given this initial condition, we can then calculate the SCB dynamics, as well as the 
outgoing field 

$™*(t) = $<_((T,t), 

propagating to the left in the line. The flux at x = is simply the sum of the incoming 
and outgoing flux fields 

$ (£) = $(Q-, t) = & n (t) + $ ou '(t) + V DC t, (10) 
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where we for simplicity also explicitly extracted the dc voltage bias Vdc-, implying that 
$ m and <& out have no dc components. Now, solving for p from (j3J) gives 

po = ^ [v DC + d t (** + «&«")] - ^ PJ , (ii) 

and inserting this expression into (jSJ), we arrive at 

= Pj + C c [V DC + d t (¥ n + ^)] ^ 

We then insert $o from ffTUl) in (TI51) and arrive at 

9iPo = z = = , (13) 

where we used the relation d x &<±(0~) = =Ff _1 (9 t $^(0 _ ), to change the spatial derivative 
into a time derivative. Inserting the expression for po from (II ip into the left hand side 
of this equation and integrating once with respect to time leads to 

^out = $i „ + Zq C^ Pj _ TRcQt (gin + Qout\ ^ (14) 

where the time trc = C c CjZq/Cy, is the characteristic RC-time for discharging the SCB 
through the transmission line. Equations (jSJ), (1121) and (1141) in principle give the full 
time evolution of the SCB operators <3>j and pj as well as the out-field, in terms of the 
in-field. However, to solve these nonlinear equations straightforwardly, we need to make 
some approximations. 

2.3. Voltage biased SCB approximation 

In the following, we will neglect the last term in ( fl4l) . Since the time-derivative enters in 
product with r RC , this will be a good approximation as long as the relevant frequencies 
of the incoming field $ m and of the SCB dynamics (pj) is much lower than the inverse 
RC-time. Under this approximation, the final equations of motion are: 



Pj + C c (V DC + 2d t ¥ n + ^dtpj^j 



dt*j = * ^ : (15) 

1e fie \ 

d tPj = - Ej- sin I -$A (16) 

^out = $ in + TlE pj . (17) 
w 

Here, we also note that this approximation is valid in recent experiments [T8l[T9] , where 
Z = 50Q, C c ~ 10 fF and Cj ~ 25 fF, giving an inverse RC-timescale of 1/(2tttrc) ~ 
400 GHz, which is around 50 times higher than the relevant frequency of $ m and pj, 
set by the qubit frequency ~ 7.5 GHz . 

The above set of equations (I15til7p correspond to the Hamiltonian 

H = H sys + H int + H bath , (18) 

„ [pj + C c V DC f w fie, \ nm 
H sys = ^ Ej cos I — $j 1, (19) 
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[c c dMo-,t)] 2 f° P (x,ty [d^(x,t)} 2 

Hbath ~ — 2ct — + -icr + 2L (21) 

Thus we have arrived at the Hamiltonian of a voltage biased SCB, weakly coupled to 
the transmission line voltage at x — 0, i.e. V (t) = d t §o(t). (Here, we note that for the 
uncoupled transmission line, without SCB, &o(t) = 2$ m (t) due to the perfect reflection.) 
From this point we can proceed with a Bloch-Redfield derivation of a master equation 
for the SCB only [34J. By comparing to section 3.2. in [30], we also note that the 
equations of motion (I15H16P can be interpreted as quantum Langevin equations (QLE) 
of the form 

Y = l j_[H sys ,Y] + ^aX - 2^A™, [X,Y}} + , (22) 
whereas (JT7J) stands for the input-output relation 

A out (t) = A m (t) - yj^X{t), (23) 

using the identifications Y = $j, X = —(pj + C c Vdc), A m = \^C^^ m and where 



-> = z °{rJ (24) 

is the damping constant that accounts for spontaneous emission. 
2.4- Master equation 

From the QLE (122]) we can derive a master equation for the reduced density matrix of 
the SCB in the transmon regime. As input field, we first consider a thermal background 
at temperature T, giving rise to a photon occupation number of 

n u = ) -. (25) 

exp [noj / k B ± ) — 1 

We assume that the density matrix initially can be written as a direct product, as well 
as Markovian properties and short correlation times for the transmission line variables. 
In the case when the damping ff24|) is much smaller than the system eigenenergies we 
arrive, after also employing the rotating wave approximation, at the following quantum 
optical master equation, 

pit) = ~ -^[Hsys,p] 

+ [("«*» + X ) V ( X ^)P + n^ m V(X+)p] , (26) 

m 

with the Lindblad operator defined by T>(c)p = cpc^ — | (c^cp + pc^c) . Also, X has been 
decomposed into eigenoperators of H sys , 

[H sys ,X±] = ±hu m X±, u m > 0, (27) 

which is always possible as long as the eigenstates of H sys form a complete set. 
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Projecting the master equation onto the SCB eigenstates H sys \i) = 
(i G {0, 1,2,.. .}), we arrive at the following equation for the diagonal elements, 

Pii ^ I" jiPjj ^ijPii (28) 

where the relaxation (wjj = Ui — uij > ) rates are 

^ij = ^ij^ + n t , ij )\(i\X\j)\ 2 , (29) 
and the excitation (cjy < 0) rates are 

Iy = yK-ln^yl (i\X\j) | 2 . (30) 

Remembering that X = — (pj + C c Vdc7) is the charge operator, the matrix elements 
can be calculated numerically from the SCB Hamiltonian in ( TT9|) . Denoting the SCB 
charging energy Ec = e 2 /2C*s, the transmon regime is found for Ej S> Ec [29]. Here, 
the SCB spectrum approaches a linear oscillator with the junction plasma frequency 
Up = y/8EjEc/h , and the charge operator asymptotically couples only neighboring 
eigenstates [29]. We find the non-zero relaxation rates 

r 0+1) , = tt(j + l)n 2 ^- + n Up ), (31) 

and excitation rates 

r i(i+i) = *U + ^^"y^^' ^ 32 ^ 

where Rk = h/e 2 ~ 25 kf2 denotes the quantum of resistance. The off-diagonal (i 7^ j) 
elements are subject to a pure exponential decay, 

Pij = -JijPij (33) 
with dephasing rates 



V M« Mi / 



7v = r ; + n + o E r ^ + E r ^ > (34) 



equal to half the sum of all rates for transitions from state |i) and as well as the 
pure dephasing rates 

K = 2 i kj f mmi2 - (35) 

The pure dephasing rates depend on the DC voltage, through the SCB spectrum, 
according to 

e du k (n q 



\(k\X\k)\ 



4E, 



c 



dn g 



(36) 



where n g = C c Vdc/^ is the dimensionless gate charge of the SCB. In the transmon 
regime the spectrum is well approximated by 

uj k {n g ) = u k (n g = 1/4) - ^cos(27rn 5 ), (37) 
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where 



fc! V 7T V.2£c/ 

giving a maximum thermal pure dephasing rate (for n g = ±1/4) of 



-nk 

max 1 ± = k 



Zp k B TTT 3 

R K h 8 



E, 



c 



(39) 



Here, we also note that in addition to small amplitude thermal charge noise there 
can also be a slow but large amplitude charge drift. In some cases, the effect of this 
drift can be taken into account by averaging over the range of transition frequencies 
involved. In the transmon regime, for the transition from \k) to \k + 1) this is given by 

e k+l — ~ e fe+l- 

2.5. Coherent drive 

In the next chapter, we will examine the scattering of coherent signals on the transmon 
in the two-level and three-level approximations. To include a coherent drive in the 
description, we take the input field $ m (t) to consist of a classical part <&™{t) on top of 
the thermal background. Deriving the master equation for this case, it turns out that 
(BE]) is modified by adding the following time- dependent term to the system Hamiltonian, 

H d (t) = -2y^^m (40) 

2.6. Adding more transmission lines 

In this section, we generalize the above master equation by adding more semi-infinite 
transmission lines to the SCB. First, by adding one more semi-infinite line, we arrive 
at the important case of an SCB capacitively coupled to an infinite transmission line. 



The discretized circuit is shown in Figure 2(a), and the corresponding Hamiltonian is 
obtained from (CQ) by adding the transmission line terms for x > 



Ax^\2C 2L, 

From a similar analysis as above, we arrive at exactly the same master equation for the 
transmon's reduced density matrix, with the replacements 

= \ (*? + *S) , r RC = f 7 = f ) 2 , (42) 

and the output fields are obtained from 

n% = + (trc/Cj) P j. (43) 

We note that the damping constant 7 as well as the RC-time trc are both halved 
compared to the semi- infinite case, since the impedance to ground is halved to Zq/2. 
The in- field is the sum of the fields incoming from the left and right, but compared 
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(b) 

Figure 2: (a) Discretized circuit describing the interaction of a single Cooper-pair 
box with microwave photons propagating in an infinite transmission line, 
(b) Generalization of the input-output formalism to an arbitrary number 
of ports connected by an artificial atom. 



to the semi-infinite case the coupling coefficient is halved, since there is (almost) no 
reflection at x — 0. Indeed, for a more general scenario with N symmetrically coupled 
incident fields, as illustrated in Figure 2(b), the mapping would be 

2 



N 



N ^ 

n=l 



Zf) C r C, 



• Trc 



J 



Zq ( C, 



N C s 



7 N V C s 



and using the relation <3> = + (Vn), the output fields are given as 

N 



out 



$0 - K 



2 

iV 



1 



c, 



(44) 



(45) 



3. APPLICATIONS: SCATTERING BY THE TRANSMON 

5.i. Two-level dynamics 

In this section, we examine the scattering of coherent signals on the transmon in an open 
transmission line. The input field is a constant coherent signal with a single frequency 
Up, close to resonance with the first transition frequency Wio of the transmon. Thus, we 
can safely describe the transmon as a two-level system. The master equation is given 
by (1261) with a coherent drive and generalized to the case of an infinite transmission 
line (see sections 12.51 and 12. 6[) . in the special case of only one system eigenfrequency 
uiq. Moreover, we include an additional term due to pure dephasing, so that the total 
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dephasing rates are given by (I54"|) . We represent our operators by the following Pauli 
matrices (using the notation X^ = (i\ X 

Hsys = -h^a z , X ± = ±i\X 1Q \a±. (46) 

Below, we will determine reflection and transmission coefficients for coherent signals 
scattered on the transmon. In the previous section the incoming and outgoing fields were 
described in terms of the flux, since that gives a simpler description of the transmon. 
However, the voltage is a more intuitive quantity than the flux and it is also usually 
what is measured in experiments. Therefore, in this section, we will describe the inputs 
and outputs in terms of the voltage. 

We consider an incoming coherent voltage field, 

V[ n (t) = n p smu p t, (47) 

impinging on the transmon from the left. For simplicity, we set the temperature to zero 
( n iu 10 = 0). The reflected voltage field is the output to the left of the transmon. Using 
flU and (jm, we have 



Vr t (t) = - ] /^{X(t)), (48) 
where the expectation value can be written as 

(X(t)) = iu w ((X + (t)> - (X~(t))) = -u l0 \X w \(a*). (49) 
Inserting (j49j) into (HE]) yields 

Vr\t) = ^hu w T 10 Z (a x ) = y/hw w T 1Q Z Re [p i] , (50) 

where poi is a density matrix element in the transmon eigenbasis. 

To solve the master equation we perform a unitary transformation to a frame 
rotating with the driving frequency u p . In this frame, the equation becomes time- 
independent after employing the rotating-wave approximation. Solving the equation 
in the steady-state (p = 0) and transforming back to the non-rotaing frame yields the 
following expression for the desired density matrix element, 



= i y/nuj w T 10 z (A + z 7l0 ) n p eiUpt ^ 

2 fiw 10 Z j% + hw w Z A 2 + 7 10 ^ 

where A = u p — uj\o is the detuning. Now, plugging this expression into (!50|) results in 

Q sin u„t — — cos ui v t 
V out (t) - p 710 (^9) 

V L W - 2 710 + + 2 > y bZ > 

Tio rio7io Tio 

where N in = Q 2 / (2Z fiiUJi ) is the average number of incoming photons per second. 
Thus, the reflection coefficient for the negative frequency part of the field is given by 
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Figure 3: Reflectance R (red) and transmittance T (blue) for a two-level transmon 
as a function of detuning, with the average number of incoming photons 
per interaction time being iVj n /(rio/27r) = 0.01. 

with r = r 10 /27 10 . For the transmitted field, f l4U]) yields 



which directly gives us the following expression for the transmission coefficient, 



In Figure [3] we plot the reflectance R = \r\ 2 and transmittance T = |t| 2 as a function of 
the detuning, in the case of a weak input signal and no pure dephasing. For a resonant 
drive (A = 0) we see that perfect reflection is approached, in agreement with [T2ll 13( 135]. 

3.2. Three-level dynamics 

In the previous section, we showed that a low-amplitude input signal is totally reflected 
when it resonantly scatters off a transmon in the two-level approximation. In this 
section, we instead study the scattering off a transmon in the three-level approximation. 
By strongly driving the second transition, the transmon becomes transparent to 
frequencies in resonance with the first transition. This effect is due to Autler-Townes 
splitting and has been observed in recent experiments [T8] . 

We consider an incoming voltage field from the left, consisting of a probe field 
flpsmojpt close to resonance with the first transition (with detuning A p = u p — uio) and 
a control field Q c sin u c t close to resonance with the second transition (with detuning 
A c = u c — u)2i)- Figure H] shows the energy levels of the transmon in this approximation. 
In the transmon eigenbasis, the relevant operators are (with the ground state energy 




(54) 



t = 1 + r 




(55) 



w = 0) 



2 




(56) 
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Figure 4: Internal levels of the transmon in the three-level approximation. A 
strong control field drives the |1) —> |2) transition, rendering a 
transparency for the |0) — > |1) transition. 



X =iJ2\Xi(i-i)\(°t-*r), (57) 
1=1 

with af = |z) (i — 1| and err = (cr J + )''\ In the same way as in the two- level case ( (JlSl) - 
( I50|) ). we obtain the following expression for the reflected signal 

V?\t) = -^f(X(t)) = \Y, y/h^m^M), (58) 

i=l 

with of = o~l + cr^ . Thus, the reflected field consists of one part with frequencies around 
the the probe frequency u p and one part with frequencies around the control frequency 
l>j c . Since we are interested in the reflectance and transmittance properties of the probe, 
we concentrate on the corresponding part of the reflected field, 

V p ref (t) = ^hw 10 Z T w (a x 1 ) = y/frw 10 Z r 10 Re(p w ). (59) 

The master equation is given by fT2oT) for the case of two system eigenfrequencies, 
again with a coherent drive and generalized to the case of an infinite transmission 
line (see sections 12.51 and [276]) . Also, terms accounting for pure dephasing are added. 
To transform the master equation to a time-independent picture, we use the following 
unitary transformation matrix, 

/ 1 \ 

U(t) = e-^f* 

y e -i(u p +uj c )t j 

and employ the rotating-wave approximation. As before, we solve the master equation 
in the steady-state to determine p w , but we now consider two different cases. 

Firstly, by setting Q c = 0, we recover exactly the same expression for the reflected 
field as in the two-level case. Thus, with the control field turned off, we see almost full 
reflection for weak probe fields on resonance with the first transition frequency of the 
transmon. 

Secondly, we consider the case of a strong control field (Q c ^> Q p ). Solving the 
master equation and expanding p 10 to first order in (Q p /fl c ), we obtain the following 
expression, 



(60) 
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Figure 5: (a) Transmittance T p of the probe as a function of detuning for 
three different control field strengths; -/V"P n /(Tio/(27r)) = 0.01 (red), 
A^/(r 10 /(27r)) = 1 (blue) and A^/(r 10 /(2^)) = 10 (green), (b) 
Transmittance as a function of control field strength for a resonant probe 
field (A p = 0). 



2itu 2 iZ J j^- (720 - i (A c + A p )) Q p 



Hl ° 4^ 21 Z (7io-iAp)( 7 20-i(A c + A p ))+r 21 02 ' 1 ; 

Inserting f lBTj) into f )59|) we can determine the reflection coefficient. For a resonant control 
field (A c = 0), the result is 

2r 10 ( 72 2 + Aj) (tip - ?Ap) + r 10 r 21 ( 720 + tA p ) N? n 
4 ( 7l 2 + a 2 ) ( 7 2 + a 2 ) + 4r 21 ( 7lo720 - a 2 ) jvf + r 2 ^ 2 1 j 

where Nf n = tt 2 c /(2Z hw 2 i) is the average number of incoming photons per second in the 
control field. The transmission coefficient is again given by t = 1 + r. Figure |5] shows 
the transmittance T = |t| 2 for different probe detunings and control field strengths. In 
these plots, we have neglected pure dephasing and used ( |29l) to express T 2 i in terms of 

r 10 . 

We can clearly see that, for strong control fields, the transmittance of a resonant 
probe approaches unity. Thus, by turning on and off a strong resonant control field, 
we can switch between the cases of full transmission and full reflection for the resonant 
probe. 



3.3. Second-order correlations 

The first-order correlation functions contain information about the intensity of the field 
and can, as we have already seen, be used to compute reflectance and transmittance 
properties in scattering processes. These functions, however, do not tell us everything 
about the photon statistics. The normalized second-order correlation function is in the 
steady state given by [36] 

(2) M _ (V + (t)V+(t + T)V-(t + T)V-(t)) 

9 {T) ~ (V+(t)V-(t)} 2 ' {b6) 
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and is proportional to the conditional probability of detecting a photon at time t + r, 
given that one was detected at time t. Here, V ± {t) are the positive and negative 
frequency parts of the voltage field. 

In this section, we calculate g^(r) for the transmitted and reflected fields from a 
transmon driven by a resonant coherent signal. We treat the transmon as a two-level 
system and use the same notation as in section 13.11 To be able to compare with recent 
experiments in [19] we perform the calculations for finite temperatures and a finite 
detection bandwidth on the output signal. For zero temperature and infinite bandwidth 
we recover the results of [37]; perfect antibunching in the reflected field and bunching 
in the transmitted field. 

Including the effect of a finite detection bandwidth is straighforward, by including 
a filter in the calculations. The approach we have taken is to model the filter by a single- 
mode transmission line resonator in resonance with the transmon, with the Hamiltonian 

H res = fyjji a)a, (64) 



and cascade it with the transmon. We start from the quantum Langevin equation ( 1221 ) 
for the transmon, generalized to the case of an infinite transmission line (see section 
12. 6p . and a similar equation for the resonator. Our coherent input signal is the voltage 
field Vl n (t) = fldsmojdt, just like in section 13.11 We can then use the formalism of 
cascaded quantum systems in [30J to arrive at a master equation for the joint density 
matrix of the transmon and the resonator. For the field reflected from the transmon, 
the result is 

i 

p = -[p, H sys + H res ] + T 10 T>(a~)p + T 01 V(a + )p 

+ iBW 



+ l)V{a)p + ^fV{<J)p 



+ ^V r ioK;io + 1)ibw ([a,po~ + ] + [a\a p)) 



V 2 (^io + 1) 

where we have denoted the filter bandwidth by ^bw- For the field transmitted through 
the transmon, f[6"o"j) is modified by simply adding the following term to the right-hand 
side, 



[p,a J -a\. (66) 

The output we are interested in is the voltage field leaking out at the right side of the 
resonator, whose positive and negative frequency parts are proportional to a{t) and a\i) 
respectively. Thus, g^(r) can be calculated as 

(2) (ot(t)ot(t + r)o(t + r)o(t)) Tr [ a t a P(r)(op B flt)] 

9 (aHtMt)) 2 Tr[ata Ps ] 2 ' 1 J 
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Figure 6: g^ 2 \r) for the field reflected from a transmon for different temperatures 
and detection bandwidths; blue (T = 0, BW= 1 GHz), green (T = 0, 
BW= 55 MHz), red (T = 50 mK, BW= 55 MHz). Typical parameter 
values from recent experiments are used (rio/2-7r = 41 MHz, ujiq/2tt = 
5.12 GHz, P=-131 dBm). 



where p s is the steady-state density matrix and P(r) the propagator super-operator, 
defined by p(t + r) = P{r)p(t). Both p s and P(r) are obtained by solving the master 
equation ( I65l) -( l66l) . For the case without filter, g^(r) for the reflected field is given by 
(157)1 with a replaced by cr~. Since (<r~ ") 2 = 0, it directly follows that # (2) (0) = 0, i.e. 
perfect antibunching. 

In Figure [6] we plot g^ 2 \r) for the reflected field for different temperatures and 
detection bandwidths. Typical parameter values from recent experiments [19] are used. 
For zero temperature and large bandwidth we see perfect antibunching, as expected. For 
a decreased bandwidth the full time dynamics of the antibunching cannot be resolved, 
which results in a less pronounced antibunching dip. For finite temperatures we see even 
less antibunching, due to a nonzero probability of detecting bunched thermal photons. 
These results explain the qualitative features of the experimental data in [19] well. 

4. Summary and conclusions 

Summing up, we have performed a thorough analysis of the qubit-photon scattering 
in a one dimensional continuum from a microscopic point of view. In particular, we 
have derived a master equation description using a superconducting transmon qubit as 
our scatterer. When we consider the two lowest levels of the transmon it behaves as a 
mirror for the incoming photons. Then, going beyond to the two-level approximation, 
we can use a control field resonant with a second transition of the transmon to suppress 
this reflection of photons at the probe frequency. Finally we discussed how the photon 
antibunching observed in the reflected field is reduced by finite temperature and finite 
detection bandwidth. 
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